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We consider Yang-Mills theories in a recently proposed family of nonlinear covariant gauges that 
consistently deals with the issue of Gribov ambiguities. Such gauges provide a generalization of the 
Curci-Ferrari-Delbourgo-Jarvis gauges which can be formulated as an extremization procedure and 
might be implemented in numerical calculations. This would allow for nonperturbative studies of 
Yang-Mills correlators in a broad class of covariant gauges continuously connected to the well-studied 
Landau gauge. We compute the ghost and gluon propagators in the continuum formulation at one- 
loop order in perturbation theory and we study their momentum dependence down to the deep 
infrared regime, with and without renormalization-group improvement. In particular, we show that 
the theory admits infrared-safe renormalization-group trajectories with no Landau pole. Both the 
gluon and the ghost behave as massive helds at low energy, and the gluon propagator is transverse 
even away from the Landau gauge limit. We compare our results to those obtained in the usual 
Curci-Ferrari model, which allows us to pinpoint the specihc effects arising from our treatment of 
Gribov copies. 
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I. INTRODUCTION 

Describing the long-distance dynamics of non-Abelian 
gauge fields is a problem of notorious difficulty, mainly 
due to the breakdown of standard perturbative tools at 
low energies. Various complementary approaches have 
been developed to tackle the problem, either through di¬ 
rect Monte Carlo simulations, or by means of approxi¬ 
mate continuum methods. The former, which use a Eu¬ 
clidean lattice discretization, can directly access physical 
observables but are essentially limited to static quanti¬ 
ties. In contrast, continuum approaches, such as trun¬ 
cations of Dyson-Schwinger equations BE], or the non¬ 
perturbative renormalization group [3], can access real¬ 
time quantities HE]. However, they essentially rely on 
computing the basic correlation functions of Yang-Mills 
fields and require a gauge-fixing procedure. It is thus 
important to have a quantitative understanding of such 
correlators in order to assess the (necessary) approxima¬ 
tions underlying such approaches. This can be achieved 
through gauge-fixed lattice calculations. 

An important issue concerns the algorithmic complex¬ 
ity of fixing a gauge numerically, in particular, in the case 
of covariant gauges, the most convenient ones for contin¬ 
uum calculations. Generically, this requires one to find 
the roots of a large set of coupled nonlinear equations, 
which is numerically demanding. However in some cases, 
the gauge condition can be formulated as an extremiza¬ 
tion problem and efficient minimization algorithms can 
be employed. This strategy was successfully applied for 


the Landau gauge leHUj. Attempts to formulate general 
linear covariant gauges on the lattice have been made in 
Refs. pMUTfi] , However, the extremization functional pro¬ 
posed in Refs. [13 El has spurious solutions, while the 
procedure of Ref. m is limited to infinitesimal deviations 
from the Landau gauge. The numerical procedure put 
forward in Ref. m and further developed in Refs. [ni- 
l2Q] seems to correctly produce linear gauge-fixed config¬ 
urations although its practical implementation does not 
really correspond to solving an extremization problem. 
Incidentally, it is problematic in what concerns the ghost 
sector. Linear gauges have been recently investigated in 
the context of Dyson-Schwinger equations in Refs. 

Ej. 

In a recent work OEll, we have proposed a simple 
generalization of the Landau gauge extremization func¬ 
tional, which generates a one-parameter family of nonlin¬ 
ear covariant gauges with a well-defined continuum for¬ 
mulation. Neglecting the issue of Gribov copies and ap¬ 
plying the standard Faddeev-Popov (FP) procedure, this 
reduces to the Curci-Ferrari-Delbourgo-Jarvis (CFDJ) 
Lagrangian [26l [27] . The latter is unitary and renormal- 
izable in four dimensions. The proposed extremization 
functional satisfies the required properties for powerful 
numerical minimization algorithms, such as, e.g., the Los 
Alamos algorithm [28]. A numerical implementation of 
this proposal would open the way for nonperturbative 
calculations of Yang-Mills correlators in a broad class of 
covariant gauges continuously connected to the Landau 
gauge. 
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The purpose of the present work is to study the ghost 
and gluon propagators in the continuum Euclidean for¬ 
mulation of this class of gauges. Continuum approaches 
have to deal with the issue of Gribov copies [29]. In the 
present context, this originates from the fact that the ex- 
tremization functional possesses many extrema. Numer¬ 
ical minimization algorithms allow one to isolate a single 
extremum (minimum), but such a procedure is difficult 
to formulate in terms of a local action, the necessary 
starting point for continuum techniques. The standard 
FP procedure completely ignores the Gribov issue and 
is valid, at best, in the high-energy regime, where one 
expects Gribov copies to be irrelevant. It is thus of key 
importance to devise gauge-fixing procedures which deal 
with the issue of Gribov ambiguities. A famous example 
is the (refined) Gribov-Zwanziger proposal [291132] . where 
one restricts the path integral to the first Gribov region, 
in which the FP operator is positive definite. This has 
been mostly worked out in the case of the Landau gauge 
but some studies exist for linear covariant gauges |33ll37] . 
This procedure is, however, not completely satisfactory 
since the first Gribov region itself is not free of Gribov 
ambiguities [38ll39] . 

An alternative strategy, put forward in Ref. [40] in the 
case of the Landau gauge, is to average over the Gri¬ 
bov copies in such a way as to lift their degeneracy. 
For an appropriate choice of the averaging procedure, 
this can be formulated in terms of a local renormaliz- 
able field theory in four dimensions, which turns out to 
be perturbatively equivalent to the Landau limit of the 
Gurci-Ferrari (CF) model [26]— a simple massive exten¬ 
sion of the FP Lagrangian—for the calculation of ghost 
and gluon correlators. Here, the bare gluon mass is re¬ 
lated to the gauge-fixing parameter which lifts the degen¬ 
eracy between Gribov copies. Remarkably, this theory 
admits infrared-safe renormalization-group (RG) trajec¬ 
tories, allowing the use of perturbation theory down to 
arbitrarily low momenta El- One-loop calculations of 
ghost and gluon correlators in the vacuum and at finite 
temperature are in good agreement with lattice results 
[ITlfe] . This has been extended to QCD in Ref. [44] 
and to background field methods in Ref. [45H47| , allowing 
for a simple perturbative description of the confinement- 
deconfinement transition in pure S\J{N) theories. 

We have generalized such an averaging procedure over 
Gribov copies to the family of nonlinear covariant gauges 
mentioned above in Ref. [24] ; see also Ref. [25] for a short 
overview. Again, for a suitable choice of the averaging 
procedure, this can be formulated as a local field theory, 
which is perturbatively renormalizable in four dimensions 
and, thus, well suited for continuum approaches. As 
in the case of the Landau gauge, the resulting gauge- 
fixed action consists of a massive extension of the GFDJ 
action—the general GF action [26] —coupled to a set of 
replicated supersymmetric nonlinear sigma models. The 
superfield sector is directly related to the treatment of 
Gribov ambiguities. However, unlike in the special case 
of the Landau gauge, this superfield sector does not de¬ 


couple in the calculation of ghost and gluon correlators 
and the resulting effective gauge-fixed theory is very dif¬ 
ferent from the general GF model. In this article, we 
present a complete calculation of the ghost and gluon 
propagators at one-loop order using this gauge-fixed ac¬ 
tion, with and without RG improvement, and we investi¬ 
gate possible infrared-safe RG trajectories. We compare 
our results to those in the GF model, which gives direct 
information on the role of Gribov ambiguities. 

One aim of the present work is to motivate lattice cal¬ 
culations of ghost and gluon correlators in the proposed 
family of gauges. Such calculations would typically be 
performed by selecting one minimum (maximum) of the 
extremization functional. It is by no means clear how this 
would be related with the averaging over Gribov copies 
studied here. However, we may hope that the situation 
is similar to the case of the Landau gauge, or at least 
close to the latter. As mentioned above, in this case, the 
weighted average procedure put forward in Ref. [40] pro¬ 
duces perturbative results which agree well with lattice 
calculations in the minimal Landau gauge—where one se¬ 
lects a unique (random) Gribov copy in the first Gribov 
region, that is, among the ensemble of minima (maxima) 
of the extremization functional. A possible interpretation 
of this observation is that, for some range of the averag¬ 
ing parameter, all Gribov copies in the first region are 
essentially equiprobable while other Gribov regions are 
suppressed, such that averaging over copies is equivalent 
to randomly choosing a particular one in the first region 
m- Understanding this issue further would require a de¬ 
tailed investigation of the landscape of the extremization 
functional. This can be done on small lattices [48] but it 
is a formidable task in general both because the number 
of copies typically grows exponentially with increasing 
volume and because it is virtually impossible to numer¬ 
ically find copies that are not in the first Gribov region 
(i.e., saddle points of the extremization functional). 

The paper is organized as follows. In Sec. [n| we recall 
the main steps of the gauge-fixing procedure proposed 
in Ref. [24]. We fix the notations and we introduce a 
crucial technical aspect of the whole procedure, namely 
the replica trick. In Sec. uni we review the Feynman 
rules relevant for one-loop computations and calculate 
the self-energies of the theory. We show in an explicit 
example how the calculation of loop diagrams involving 
the superfield sector are performed in practice. In Sec. 
|TV| we present two renormalization schemes and we dis¬ 
cuss, in particular, the role of the replica. In Sec. [V| we 
present our results for the ghost and gluon propagators 
for different choices of the gauge-fixing parameter. Fi¬ 
nally, we investigate RG improvement in Sec. |VI| and we 
find infrared-safe RG trajectories where the coupling con¬ 
stant remains under control all along the flow. This justi¬ 
fies the use of perturbation theory down to deep infrared 
momenta. In this regime, we find that both the gluon 
and the ghost are massive, and that the gluon propaga¬ 
tor is always strictly transverse. The superfield sector is 
crucial in that it ensures that the present approach is in- 
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deed a bona fide gauge fixing, in contrast, e.g., to the CF 
model. In order to emphasize the role of the superfield 
sector, we systematically compare our results to those of 
the CF model at the same order of approximation. Fi¬ 
nally, we give some technical details in Appendices [A[|D| 


II. THE GAUGE-FIXING PROCEDURE 

The classical action of the S\J{N) Yang-Mills theory 
reads, in d-dimensional Euclidean space, 

= 4 y ’ (1) 


where f^ = f and 

f ;, = - d,Al + gor'^’^AlAl, ( 2 ) 

where go is the (bare) coupling constant and a summation 
over repeated spacetime and color indices is understood. 
In the following, we use the convention that fields written 
without an explicit color index are contracted with the 
generators of SU(Y) in the fundamental representation 
and are thus N x N matrix fields, e.g., = A^t^. Our 

normalization for the generators is such that 


Ft'’ 



j^fabc ^ ^abc 


( 3 ) 


where and are the usual totally antisymmetric 
and totally symmetric tensors of S\J{N). In particular, 
we have 


tr (Ft'’) 


gab 

~Y' 


( 4 ) 


In order to fix the gauge, we consider the functional 


niA,g,U]= [ tr 

J X 


jj^2 ^ U'^rj + ig^U 


«) + 


( 5 ) 


where 77 is an arbitrary N x N matrix field, and 

= UA^U^ + 


(6) 


is the gauge transform of A^ with U G SU(Y). We define 
our gauge condition as (one of) the extrema of 1-L with 
respect to U. This leads to the following covariant gauge 
condition: 


(5^0“ = ^tr[F([/7?t-r,[/t)]. (7) 

The functional H admits many extrema Ui = Ui[A^ri\ 
which correspond to Gribov copies. The possible nu¬ 
merical implementation of this gauge condition has been 
discussed in Ref. m- Once the gauge condition 0 has 


been implemented, we average over the random field 77 , 
with a Gaussian weight^ {M is a normalization factor) 


V [rj] = AT exp 



( 8 ) 


where ^0 is a (bare) gauge-fixing parameter. The case 
^0 = 0 corresponds to 77 = 0, that is to the Landau gauge 

d,A\i = 0 . 

Let us briefly recall the main lines of the formalism 
developed in Ref. [24]. We define the vacuum expecta¬ 
tion values of an operator 0[A] by a two-step averaging 
procedure. The first step, hereafter denoted with brack¬ 
ets, consists in an average over the Gribov copies of any 
gauge held configuration A: 


' ^ " JVgV[7i]J2iS(i)e-0o-H[A,^,u,] ’ 


( 9 ) 


where the discrete sums run over all Gribov copies, s{i) 
is the sign of the functional determinant of the Faddeev- 
Popov operator—the Hessian of the functional 0 
evaluated at U = Ui^ and /3o is a free gauge-fixing pa¬ 
rameter which controls the lifting of degeneracy between 
Gribov copies. The case /3o = 0 would correspond to a 
flat weight over the Gribov copies and would reduce to 
the standard FP construction. The second step, denoted 
by an overall bar, is a standard average over gauge held 
configurations with the Yang-Mills weight. The expecta¬ 
tion value of the operator 0[A] is thus obtained as 


{0[A]) 


/F>Ae--SYM[^] 


( 10 ) 


Using standard techniques, the sum over Gribov copies 
in Eq. ^ can be written in terms of a path integral over 
an SU(Y) matrix held U as well as ghost and antighost 
fields c and c and a Nakanishi-Lautrup held h which en¬ 
sure the gauge condition 0 . One can then explicitly 
perform the Gaussian integration over the held 77 . A 
crucial point here is the presence of the denominator in 
Eq. 0 , which guarantees that gauge-invariant observ¬ 
ables Oinv[^ are blind to the gauge-fixing procedure: 
{Oinv[M) = Oinv[^‘ Such a denominator produces a 
nonlocal functional of the gauge held which can, however, 
be treated by means of a local action using the replica 
technique m- The latter amounts to making n replicas 
of the set of fields (U, c, c, h) and taking the limit 77 ^ 0 
at the end of (loop) calculations. One can then explicitly 
factor out the volume of the gauge group f VU, which 
effectively singles out one replica. The 77 — 1 remaining 
replicated fields (U/c, c/c, c/^, h/^), with 2 < /c < 77, are con¬ 
veniently grabbed together by introducing SU (N) matrix 


^ The weight ^ is invariant under the change 77 —)■ r]U, which 
ensures that one can eventually factor out the volume of the 
gauge group in the FP construction. 
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superfields Vk that depend on the Euclidean coordinate 
X and a pair of Grassmannian coordinates {Ok^Ok) = 
as 

= exp ^igo{9kCk + CkOk + (H) 

where hk = ihk + c/c}. Here, the x dependence on 

the right-hand side is only through the fields. Note that 
there is an independent pair of Grassmannian coordinates 
for each replica. 

After some manipulations, one finally gets m 
f V{A, c, c, h, {V}) 0[A] 


{0[A]) ^im jx>(A,c,c,h,{V})e-siAc,c,h,{v}] ’ 

( 12 ) 

where P(A, c, c, h, {V}) = V{A, c, c, h) x Y[k= 2 ^^k and 
the gauge-fixed action for arbitrary n is given by 

S[A, c, c, h, {V}] = 5 'ym[^] + c, c, h] 

n 

+ ^ 5'sUSy[^, V/e], (13) 

k=2 

where 5 'ym[A] is the Yang-Mills action. 


Scy[A,c,cM = J 

{h ) go nabc’-ua-b^c 9o 


+ ^0 
+ /^O 


2 2 4 ^ ’ 




\ }■ 


( 14 ) 


and 


5'susy[^, V/c] = Scf[A^^ ,Ck,Ck,hk]. (15) 


It is remarkable that the integration over the random ma¬ 
trix field 77 with the weight ^ produces the GFDJ ac¬ 
tion, the first two lines on the right-hand side of Eq. (14). 


The third line, which corresponds to a bare mass term 
for both the gluon and the ghost fields proportional to 
the factor ^do, arises from the average over Gribov copies. 
Altogether, the first line of Eq. (13) corresponds to the 


general CE action, whereas the second line is the contri¬ 
bution from the replicated superfields. 

Einally, for actual calculations, it is useful to exploit 
the superfield formulation Eq. (© and to rewrite ( p!^ 


as 


5susy[AV] = 4 / tT{D^V^D^V+^g^^dNV^dMV}, 

9o JX^o_ 

(16) 

with D^V = d^V^-igoVA^ and where the uppercase latin 
letters M, N stand for Grassmann variables 0,0 {a sum 
over repeated indices is understood). Here, we have in¬ 
troduced a curved Grassmann manifold with line element 
ds^ = gMNdNdM = 2gQQd0d0, where 


9ee — ~9ee — PodO -h 1 , 

gOO = _g00 = _ 1 . 


(17) 


The curvature of the Grassmannian manifold is con¬ 
trolled by the lifting parameter f3o- As already noticed, 
there is such a curved Grassmannian manifold associated 
to each replica superfield Vk- The integration measure is 
defined accordingly as m 


where 


J^ = I d9d9g^/\e,e), 


g^/‘^{9,9)=g^/^{9)=M9-l. 


(18) 


(19) 

The formulation ([T^ makes transparent a large class of 
supersymmetries which mix the original bosonic (Uk^hk) 
and Grassmannian {ck,Ck) fields. These simply corre¬ 
spond to the isometries of the Grassmannian manifold 

El [50]. 

We stress again that the limit n —> 0 in Eq. ( 12 ) results 
from the presence of the denominator in Eq. li and is 
thus crucial in order for the present theory to correspond 
to a gauge-fixing procedure, as discussed above. Instead, 
ignoring the denominator in Eq. ^ is equivalent to set¬ 
ting n = 1 which does not correspond to a gauge-fixed 
theory. In that case, there is no superfield sector and the 
model reduces to the GE action, given by the first line 
of Eq. (13). It is, therefore, interesting to compare the 
results obtained in the limit n ^ 0 to the corresponding 
ones in the CE model, obtained by setting n = 1 in order 
to highlight the importance of the superfield sector and of 
the limit n ^ 0. This allows us to study the influence of 
the Gribov ambiguities. It is worth emphasizing that the 
action (13) is perturbatively renormalizable and asymp¬ 
totically free in four dimensions for an arbitrary value of 
n despite the presence of the replicated matrix superfield 
El • We can thus use the same renormalization schemes 
in the cases n ^ 0 and n = 1 . 

To summarize, we have presented a genuine gauge¬ 
fixing procedure of the Yang-Mills action in a certain 
class of nonlinear covariant gauges, which consistently 
deals with the issue of Gribov ambiguities and which can 
be written as a local, renormalizable gauge-fixed action. 
These gauges are indexed by a parameter and the 
lifting over the Gribov copies is controlled by an extra 
gauge-fixing parameter /3o > 0 homogeneous to a square 
mass. The resulting gauge-fixed action (13) depends on 


the gauge field A, the usual ghost and antighost fields 
c, c and the Nakanishi-Lautrup field h, and a set of n — 1 
replicated SU(Y) matrix superfields. In the following, we 
compute the two-point vertex functions of the theory in 
the vacuum with the action (13) for fixed n, from which 
we obtain the various propagators after inversion. We 
then either perform the limit^ n ^ 0 as required by our 
gauge-fixing procedure, or set n = 1 in order to compare 
with the CE model. 


^ This has to be done after inversion of the matrix of two-point 
vertex functions because the latter involves summations over the 
replica. 
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III. PERTURBATION THEORY 

We compute the ghost and gluon propagators at one- 
loop order in perturbation theory. The calculations are 
lengthy but straightforward. Here, we solely display the 
relevant Feynman diagrams and detail one explicit ex¬ 
ample involving the superfield sector. The complete ex¬ 
pressions of the diagrams are given in the supplemental 
material m- 


where S{0^0') = {0 — 0'){0 — O') is the covari¬ 

ant Dirac delta function on the curved Grassmann space: 
Jq ^(^5 ^0/(0 = Notice that, for 7^ 0? there is a 

nontrivial correlation between different replica. Finally, 
there are nontrivial mixed correlators 

^ (26) 

and 


A. Feynman rules 


[Ki-p^^)K(p)]o 


p^{p^ + MoY 


(27) 


For the purpose of the perturbative expansion, we in¬ 
troduce the replicated superfields A^{x,0j^) as 


Vk{x,0k) = exp [igo^Alix,^;^)]. 


( 20 ) 


Expanding the action I® in powers of up to order 
^ 0 , we get the free propagators and the vertices relevant 
for one-loop computations. We work in momentum Eu¬ 
clidean space with the Eourier convention —iPfi- 

Since the Grassmann spaces are curved, it is not useful 
to introduce associated Eourier variables. 


The quadratic part of the action (13) couples the dif¬ 
ferent fields, which leads to mixed propagators. In the 
(A, c, c, h) sector, this yields the following expressions for 
the bare tree-level propagators [24] (see also Appendix 


0: 


[A“(-p)4(p)], = V'' 


p^uip) CoP^:u{p) 
+ nYo + /3oCo 


( 21 ) 


The interaction vertices are obtained from terms higher 
than quadratic in the fields. Erom Eq. ( [T^ , it appears 
clearly that the vertices of the sector (A, c, c, h) are iden¬ 
tical to those of the CE model or, equivalently, to those of 
the CEDJ action. These include the Yang-Mills vertices 
with three and four gluons as well as the standard gluon- 
ghost-antighost vertex. In addition, there is a four-ghost 
vertex as well as a hcc vertex, both proportional to ^o- 
These vertices are well known [52] and we do not recall 
their expressions here. The vertices of the replicated non¬ 
linear sigma model sector are obtained by expanding the 
exponential (20) in powers of A/.. We refer the reader to 
Ref. [24] for details. Eor illustration, we simply recall the 
expression of one cubic vertex involving two superfields 
A\{p 2 , 0 ^) and one gauge field A^^{ps): 

i^f°-^<^Ski{2TTYS^'^Ypi+P2+P3)S{0,6'){pi -P2)m- (28) 


and 


£ab 

(22) 


(23) 


(24) 


Here, P^^(p) -Pjj/(p) = ^k-i^-PfiPu/p'^, and the 

square brackets with subscript 0 denote an average with 
the quadratic part of the action ( [T^ , with n finite. The 
only reminiscence of the replica sector of the theory in 
the correlators (21)-(24) is through the square mass nf3o 
appearing in the transverse part of the gluon propagator. 
Eor n = 1 we recover the tree-level propagators of the GE 
model. 

The correlator of the superfields Ak reads 


[Ali-p,e)A’i{p,^)l 


= (5“'’ 


Ski5{e,0') 
p‘^ + 


Co(l + 5kl) 


p^ijp^ + Mo) \ 

(25) 


B. Two-point vertex functions at one-loop order 

We are primarily interested in the ghost and gluon 
propagators. The latter are obtained from the two-point 
vertex functions, defined as the second derivatives of the 
effective action F at fixed n\ 




6X{p,e)6Y{-p,g) ’ 


(29) 


where X and Y denote any of the (super)fields, the sub¬ 
script 0 means that the derivative is evaluated at van¬ 
ishing fields, and Sq/SX is a covariant functional deriva¬ 
tive m SO] which takes into account the curvature of 
the superspaces associated with each superfield. It is 
defined as 5q/5X{p) = S/SX{p) for normal fields and 
6o/6X{p^O) = /SX{p,0) for superfields. In par¬ 

ticular, one has SoX{0)/6X{^) = 6{0,0') for a given su¬ 
perfield X. 

All two-point vertex and correlation functions are di¬ 
agonal in color space and we extract a trivial unit matrix 
from their expressions below. Spacetime symmetries im¬ 
ply the following general decompositions: 

Ta.aAp) = PUp)^t{p) + P^M^Lip) (30) 
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and 


^ihA^ip) = -^A^ihip) = ip^.I'ihAip), (31) 

where the scalar functions Tt{p), Tl{p), and TihA{p) 
only depend on p‘^. For the superfield sector, the replica 
permutation symmetry and the isometries of the Grass- 
mann subspaces associated to each replica imply the fol¬ 
lowing general decompositions: 


rA,A,{p,e,0') = dki [ri{p)6{e,0:) + r2{p)aej{e,g)] 

+ {6ki - l)rs{p) (32) 

and 

^A^A^{p,d) =-ipAiip), (33) 

where the scalar functions ri_,,.^ 4 (p) only depend on 
and where is the Laplace operator on the curved 
Grassmann space; see Appendix One has the identity 
= —2 + ‘2/3o6lo,0'). The various components 
of the two-point vertex functions at tree level are given 
in Appendix [B) 


The calculation of propagators from the two-point ver¬ 
tex functions is complicated by the fact that the fields , 
ih and A/, are mixed at the quadratic level in the effective 
action. It is, in particular, the reason why we must take 
the limit n ^ 0 only after having inverted the matrix 
of vertex functions in field space. This is discussed in 
Appendix!^ Because this mixing only involves the lon¬ 
gitudinal component (with respect to momentum) of the 
gluon field, the ghost and transverse gluon propagator 
are given by the simple relations 

G'gh(p) = lim ry(p) (34) 


and 


Grip) = lim r^^(p), (35) 


where we have decomposed the gluon propagator (p) 
as in Eq. (30). In contrast, the longitudinal component 
of the gluon propagator involves various vertex functions 
in the sector A/.) (see Appendix [A|): 


Gl{p) = lim 


^ihih (hi + /3o (n — 2) Ts) 


n^o {TLTihih — {TihAyp‘^) (hi + /3o (n - 2) Ts) - TihihP‘^l^o(^4)‘^{n — 1) 


(36) 


Note that in the case n = 1, one recovers the stan¬ 
dard expression of the longitudinal propagator in the GF 
model 


gCF 


(p) = 


- ihih 


(p) 


rL{p)rihih{p)-P^ [rihA{p)Y 


(37) 


n=l 


Using Slavnov-Taylor identities for the (non-nilpotent) 
Becchi-Rouet-Stora-Tyutin (BRST) symmetry of the GF 
model (see Appendix [C|) , this can be rewritten as 


GTip) = - 


ihih (p) 


rL(p)rcc(p) 


(38) 


n=l 


We display in Figs. the relevant Feynman dia¬ 
grams for the one-loop calculations. We use the stan¬ 
dard graphical conventions for the gluon (wiggly) and 
ghost (dashed) lines. The plain line represents the super- 
field correlator (25). The second diagram on the second 


line involves the mixed A-A correlator (27). Double lines 
stand for the field ih. 

The calculation is lengthy but straightforward. In or¬ 
der to illustrate a typical example involving the Grass¬ 
mann subspaces, we treat explicitly the first diagram 
on the second line of Fig. which involves the vertex 
Eq. (28) and a loop of superfields. According to the 


Feynman rules, its contribution to the gluon self energy 
is given by 



FIG. 1: One-loop diagrams for the vertex Vaa- The diagrams 
of the first line are present in the Landau gauge and a fortiori 
in the CF model. The diagrams of the second line involve the 
superfield sector and are thus specific of the present gauge 
fixing (they are proportional to n — 1 and thus vanish in the 
CF model). The first one is proportional to /3o^o and the 
second one to /3o^o- 
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-► 



►- 







► 








FIG. 2: One-loop diagrams for the vertex Fee- Only the first 
diagram on the first line is present in the Landau gauge. All 
the others are present in the CF model. Note that the dia¬ 
grams on the second line involve the mixed ih-A correlator 
There is no diagram involving the superfields. 


* s \ 





FIG. 3: One-loop diagrams for the vertices Tihih (left) and 
^ihA (right). Both diagrams are present in the GF model. 
There is no diagram involving the superfields. 


t~\A— loop 

j- 


(p) = 


Oil P P 

= -N / / 




hi=2 

'k 


X J 


(l+%)Co 


[ k^ + /?oCo 

fc2(fc2 

X J 


(1 + 5ij)Co 


[ + Mo 



\ 


where a, b and /i, u are respectively the color and Lorentz 
indices of the external gluon legs carrying momentum p 
and — p^. We use the notation f^=p^f 

with d = 4 — e, where the arbitrary scale p is introduced 
for dimensional reasons. The replica indices i^j are asso¬ 
ciated with the internal superfield lines. 

The Grassmannian integrals are trivially performed us¬ 
ing the identities 

[ S{0,0')f{0) = f{0 '), 6{0,0) = 0 and / 1 = /^q. 

Je_ Je_ 

(40) 

After summing over the replica indices, we obtain the 
following momentum integral: 

G:iT(p) = -5“"(n - 

HU Q 

f {k^{kj^iy) /2 2 (n + 2)^oCo\ 

"" Jk (k^ + Mo) (i^+Mo) W ^ J ’ 

(41) 


which is logarithmically divergent in the ultraviolet (UV). 
A similar calculation yields, for the second diagram on 



FIG. 4: One-loop diagrams for the vertex Faa. 





FIG. 5: One-loop diagrams for the vertex Faa- 


the second line of Fig.[^ also involving the replica sector. 




I 


d-1 
k^P^ — {k • p)‘^ 




(42) 


It is transverse and UV finite. As expected, both replica 
contributions are oc n — 1 and vanish identically in the 
CF model (n = 1). 

The calculation of the other diagrams of Figs. □"[llg oes 
along the same lines and we are left with one-loop mo¬ 
mentum integrals, whose explicit expressions are given in 
the supplemented material m- These integrals are UV 
divergent and can be regularized by standard techniques. 
We use dimensional regularization with d = 4 — e. It was 
shown in Ref. [24] that these divergences can be elimi¬ 
nated at all orders of perturbation theory by means of 
six independent renormalization factors, as we now re¬ 
call. 


IV. RENORMALIZATION 


We define renormalized fields in the standard way 




and 


= v^4, c“ = v^4, 

(43) 


h^ = ^hh\, A“ = v^A^. 


(44) 


In the following, we only refer to renormalized fields and 
correlators and we suppress the index R for simplicity. 
A nontrivial issue concerns the interplay between renor¬ 
malization and the limit n ^ 0. As pointed out in 
Refs, [m 00], part of the explicit n dependence of the 
bare theory can be absorbed in the definition of renormal¬ 
ized parameters. For instance, one observes from Eq. (21) 


that the tree-level square mass of the transverse gluons is 
n/3o. Hence, defining renormalized gauge-fixing param¬ 
eters /3 and ^ as /3o = Z^/3 and makes the 

(transverse) gluon massless in the limit n ^ 0. In the 
case of the Landau gauge (^o = 0)? the theory thus re¬ 
duces to the standard Faddeev-Popov Lagrangian m- 
















Instead, one may define a renormalized squared mass pa¬ 
rameter m? as npo = such that the transverse 

gluon propagator remains massive in the limit n ^ 0. 
For ^0 = O5 it was shown in Ref. m that the resulting 
theory is perturbatively equivalent to the Landau limit of 
the CF model. As already mentioned, the latter admits 
infrared-safe renormalization-group trajectories [41] and 
provides a perturbative description of Yang-Mills corre¬ 
lation functions in good agreement with lattice results, 
both in the vacuum muniii] and at finite temperature 
[43] . This motivates us to consider similar renormaliza¬ 
tion schemes in the case ^0 7 ^ 0- 

In that case, however, the parameter does not only 
appear in the transverse gluon mass, but also in the com¬ 
bination in the tree-level square masses of the longi¬ 
tudinal gluon and of the fields (ih, c, see Eqs. (21 )- 

(27). In order to keep this mass finite in the limit n ^ 0, 
we choose the following definitions of the renormalized 
mass and gauge-fixing parameters^: 

n/3o = Z^ 2 m‘^ and — = (45) 

such that . Important consequences of 

this choice are that the ghost propagator acquires a non¬ 
vanishing mass and that various correlators, e.g., which 
are proportional to vanish in the limit n ^ 0 at fixed 
For instance, it is easy to show diagrammatically that 
^ihih{p) oc ^0 up to a function of the parameters n^o and 
/^oCo- One thus has Tihih{p) ^ 0 in the limit n ^ 0 at 
n^o and /3o^o fixed. The function TihA{p) being finite in 
this limit, we conclude that the longitudinal gluon prop¬ 
agator (36) vanishes identically 


(46) 


Gl{p) = 0. 


We stress that this is a consequence of the definitions 1 
and that it is independent of the set of renormalization 
prescriptions discussed below. This is a remarkable prop¬ 
erty of the present gauge fixing since an exactly trans¬ 
verse gluon propagator is usually expected to be a pecu¬ 
liar feature of the Landau gauge. In particular, the fact 
that the gluon propagator remains transverse at ^ 7 ^ 0 
is a nontrivial effect of the average over Gribov copies. 
This is one of the main differences with the CF model, 
where the longitudinal gluon propagator is nonzero and 
is given by Eq. (37). 


Einally, we define the renormalized coupling g as 
go = Zg g. 


(47) 


^ This choice is arbitrary. One could also choose to have a di¬ 
verging mass for these fields, in which case they would decouple 
and lead to a different phenomenology. Lattice data would be 
needed to determine the most relevant scenario. We stress that 
such definitions of the renormalization constants are consistent 
with the renormalizability of the theory, since the proof of renor- 
malizability given in Ref. m holds for generic n. 


It was shown in Ref. [24] that the combination Z^ 2 Zc/Zh 
is UV finite so that only six renormalization factors, out 
of the seven introduced so far, are needed to absorb all 
the UV divergences. The finite parts are to be deter¬ 
mined by means of specific renormalization conditions. 
In the following subsections, we discuss two sets of pre¬ 
scriptions for the two-point functions and we use a gen¬ 
eralization of the Taylor scheme for the coupling. We 
emphasize that all our prescriptions are set before taking 
the limit n ^ 0 , which is done only at the end of the 
computations. 

As an illustration, let us consider the superfield contri¬ 
butions to the transverse gluon two-point vertex in the 
limit n ^ 0. We have, from Eq. (41), 


lim F 

n^O 


A-ioop/ _ 9 N^m y 

{d-l) 


{P) 


i 






1 1 

^ + ^ + (48) 


whereas the other superfield contribution (42) is propor¬ 
tional to oc and vanishes in the limit n ^ 0 . 


We stress again that the contribution (48) is absent in 


the CF model and is a direct effect of the average over 
Gribov copies. For instance, the contribution of the su¬ 
perfield sector to the transverse gluon two-point vertex 
at zero momentum is 


lim F: 

n^O 


A—loop 


(p = 0) = 


g^N^w? 

d 

g^N^w? 

327r2 


L 


2k‘^ + 




(49) 


where we neglected terms 0{e) in the last equality and 
p? = 47 re“^/i^, where 7 is the Euler constant. 

Another remarkable consequence of the prescriptions 


(45) is that despite the presence of a nonzero gauge-fixing 


parameter the theory possesses some of the properties 
of the Landau gauge. The fact, already discussed, that 
the gluon propagator is exactly transverse is one example. 
But there are other similar features in the loop contribu¬ 
tions. In particular, various one-loop diagrams which are 
(X ^0 vanish in the limit n ^ 0, just like the superfield 
contribution (42) to the gluon two-point vertex discussed 


above. Eor instance, it is easy to check that all but the 
first contributions to the ghost two-point vertex depicted 
in Eig. are proportional to /3o^o ^ and thus 

vanish in this limit. The first diagram of Eig. is al¬ 
ready present in the Landau gauge—with the difference 
that, here, the internal ghost and gluon propagators are 
massive—and is easily shown to be proportional to . 
We conclude that, in the limit n ^ 0, the ghost two- 
point vertex function at vanishing momentum does not 
receive correction at one loop, that is. 


lim Vccip = 0 ) = f lim Zc 
n^O ’ Vn^O 


ZpZrr, 


2 ^ . (50) 
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A similar argument shows that the one-loop diagrams of 
Fig.j^also vanish in the limit n ^ 0. More precisely, they 
are respectively proportional to and to ^o- If follows 


limVihAip) 

n^O 


hm 

n^O 


and 


lim ^ 

n^O Tiq n^O 


Z^Zh and Z^Z^ are UV finite in the limit n ^ 0. 

Finally, we employ the same renormalization prescrip¬ 
tions as described above for the CF model. In that case, 
n = 1 and there is no issue with the n dependence of 
the renormalized parameters. Note also that the non- 


renormalization relations (50)-(52) do not hold in that 


case. Instead a relation which is known to hold at all or¬ 
ders in perturbation theory in the CF model is that the 
combination (Z|ZcZ^ 2 /Z^)n=i is finite [5Ql l53l l54] . 


A. Zero-momentum scheme 


In this first renormalization scheme, we define the two 
renormalized masses from the transverse gluon and the 
ghost two-point vertex functions at zero momentum, that 
is, 


^t{0) = and rcc(O) = 


In order to fix the remaining renormalization factors, we 
further require that the following vertex functions assume 
their renormalized tree-level expressions at a scale /i^: 



= 

(54) 

Tcc{p) 


(55) 

ri(M) 

= P, 

(56) 

^ ihihip) 

= —n^. 

(57) 


This leads to the following expressions, defining dZ^ = 
Zr^-l: 


, i.e.. 

SZa = 

^t{p) — nT(o) 

5 

(58) 

(51) 

II 

szJ 

\ 

1 + 4 ) 

J 

(59) 


II 

Uccip) - ncc(o) 


(60) 


II 

.V/ 

ncc(O) rry 

irn? 


(61) 

(52) 

SZa = 

ni(M) 

’ 


(62) 

Zm‘^ , 

SZh = 



(63) 


where the various 

functions I1t{p) 

, ncc(p), 

etc. de- 


note the loop diagram contributions to the corresponding 
(renormalized) two-point vertices with the same index.^ 
Extracting the divergent parts Z^^^ oc 1/e and denoting 
K, = g‘^N/S7r‘^e, we get [24] 


SZf^ = 


13 


13 nA ^ 

6 2 J 6 ’ 








SZf^ = 




3 

4 4 

35 < 

-h — 

12 4 

~6 ~ T 

3 

4~ 12 
13 


T’ 


35 

~12^’ 


= 


13 

K ->■ —K, 

3 

' ^ 4^’ 


(64) 

(65) 

( 66 ) 

(67) 

( 68 ) 

(69) 


where we indicated the various n ^ 0 limits. We check 
that the nonrenormalization relations mentioned below 


gauge resu. 
have SZ^^^ 


(53) Eqs. (50)-(52) are satisfied. We also recover the Landau 


ts of Ref. [40] in this limit. In particular, we 
5Z^^^. Eor n = 1, we check that the above 
divergent parts coincide with the known ones of the CE 
model laisaiss], except of course for Z\. 

As for the finite parts of the renormalization factors, 
note that, in the case n ^ 0, the relation (50) combined 


with the second prescription in Eq. (53) implies that 


lim ZcZ^Z ^2 = 1. 

n^O 


(70) 


Similarly, the relation (52) and the prescription ( |57[ ) yield 
lim ZhZ^ = 1. (71) 


n^O 


We have checked this explicitly at one-loop order but we expect 
the general argument developed here to be valid at any loop 
order. 


For instance we have rx'(p) = + Z^ 2 nn?) -h 
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For later use, we also note that the relation (50) is equiv¬ 


alent to liuin^o ncc(O) = 0 and thus implies that 


lim Zc 

n^O 


ncc(/^) 




(72) 


n=0 


We also check from Eq. (|61| that Eq. (70) is indeed sat¬ 


isfied at one-loop order. Similarly, the relation (52) is 


equivalent to Uihih {p)/{n^) = 0 and we check, 

from the prescription (63), that Eq. ( [7T| is indeed satis¬ 
fied. 


B. Infrared-safe scheme 

Our second set of renormalization prescriptions is in¬ 
spired from the infrared-safe scheme put forward in 
Ref. [41], that we adapt to the case ^ 7 ^ 0. As we shall 
see below, this also leads to infrared-safe RG trajectories 
in the present case. In this scheme, all renormalization 
prescriptions are defined at the scale fi. We keep the set 
of prescriptions Eq. (54)-(57) of the previous scheme, and 


replace the prescriptions (53) at zero momentum by 


= 1 and ZlZcZ ^2 = 


(73) 


This leads to similar expressions as Eqs. (58)-(63), which 


we do not write explicitly. The divergent parts of the 
counterterms at one-loop order are the same as before 
[see Eqs. (64)-(69)], as required for renormalizability. As 
already mentioned, in the CE model (n = 1), the sec¬ 
ond relation in Eq. (73) is true for the divergent parts at 


all orders of perturbation theory [snusaiHi- Eor arbi¬ 
trary n, this relation is satisfied by the divergent parts 
at one-loop order, as can be checked from Eqs. (64)- 


(69). Here, we use this relation to define the finite parts 


of the renormalization factors. We stress that although 
we do not know whether this relation is compatible with 
renormalizability—that is, whether it is satisfied by the 
divergent parts—at higher orders (except for n = 1 ), it 
is always possible to use it to define the finite parts.^ 
Obviously, in this renormalization scheme, the inverse 
gluon and ghost propagators at zero momentum are not 
directly given by the parameters m? and anymore. 

Einally, we note that, in the limit n ^ 0, the rela¬ 
tion ([tT]), which follows from Eqs. (I^ andj^, remains 


valid. In the present scheme, the relation (51) combined 


with the first prescription in Eq. (73) implies that 


lim ZhZA = 1 . 

n^O 


(74) 


Combining this with Eq. ( pT] ) and with the second pre¬ 
scription in Eq. (73), we obtain that lim^^o Za/Z^ = 1 


and, hence, that the relation (70) is also valid in the 
infrared-safe scheme.^ 


The relations (51) and (52) also imply that the renor¬ 


malization factor Zc is identical in both schemes. Indeed, 
in the present scheme we have 




ncc(/i) ^ 






^ihih {p) 


Zi^ihA{p) 


(75) 


The relations ( [^ and (52) imply lim^^o (p) = 
\imn^oIlihih{p) = 0 and Eq. (72) holds here too. 
The (combination of) renormalization factors Zc and 
ZcZ^Z^ 2 —the only ones which enter the calculation of 
the ghost propagator at one-loop order—are thus iden¬ 
tical (in the limit n ^ 0 ) in the two renormalization 
schemes considered here. We conclude that the latter, as 
obtained from strict perturbation theory at one-loop or¬ 
der, is identically the same in the two schemes provided 
one uses the same value for the parameters m^, and 
g. Strictly speaking, this remark does not hold when RG 
improvement is taken into account since the running of 
the various parameters are different in the two schemes. 


C. Renormalization of the coupling constant in the 
Taylor scheme 

We follow the general strategy of Ref. [56] to renor¬ 
malize the coupling constant. We fix the value of Zg 
[Eq. ( |47| )] from the ghost-antighost-gluon vertex at van¬ 
ishing ghost momentum, r^cc(p, 0, —p). The Eeynman 
diagrams contributing to Tacc ^.t one-loop order are 
shown in Eig. m The two usual Yang-Mills diagrams (first 
line of Eig. [ 6 ) involve, at the vertex where the external 
ghost leg is attached, a Lorentz contraction between a 
gluon propagator and the (internal) antighost momen¬ 
tum. Eor vanishing external ghost momentum, both the 
antighost and the gluon propagators in the loop carry 
the same momentum (up to a sign). It follows that only 
the longitudinal gluon propagator contributes to these 
diagrams which are, thus, proportional to ^0 oc —up 
to a function of = Z^Z^ 2 ^m‘^ [see Eq. (21)]—and 
vanish in the limit n ^ 0. It is easy to check that all the 
other diagrams are also proportional to ^0 oc up to a 
function of and also vanish in this limit. We con¬ 
clude that the vertex Tacc{p^ 0 , —p) does not receive any 
radiative correction at one-loop order in the limit n ^ 0 . 
The divergent parts of the renormalization factors thus 
satisfy the relation 


lim ZaZc\/ZA = 1 - 

n^O ^ ^ 


(76) 


® In general, one can always replace the second relation in Eq. ( |73| 
by Z‘^ZcZ ^2 = ZaZ^ with Z = l-\-0{g‘^) chosen such that ^^ 2 , 
Zai Zct and Z^ have the correct divergent parts. 


We mention that the combination of the relations derived here 
yield lim^^^o ZaZcZ ^2 = 1, which corresponds to the infrared- 
safe renormalization prescription proposed in the case of the Lan¬ 
dau gauge in Refs. Esmi]. 
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FIG. 6: One-loop contributions to the ghost-antighost-gluon 
vertex function Facc- The diagrams on the first line are 
present in the standard Landau gauge—here with massive 
gluon and ghost propagators—and are proportional to the 
antighost external momentum. The topologies on the sec¬ 
ond and third lines are present in the standard CFDJ gauge 
and give contributions proportional to ^o- The replicated su¬ 
perfield sector only contributes in the last diagram, which is 
proportional to and to the antighost external momentum. 


nonrenormalization relation 


rv rv rv^ _ rv‘^j^ 


(78) 


which is valid for the divergent parts. As discussed above, 
we extend this equality to the finite parts, which gives us 
a definition of the renormalization coefficient Zg. It can 


be checked that this definition coincides with Eq. (76) in 
the limit ^ ^ 0 . 


V. ONE-LOOP RESULTS 

We present om results for the gluon and the ghost 
propagators at one-loop order in the SU(3) theory^ 
within the two renormalization schemes discussed above. 
In the following, we set the scale fi = 1 GeV and we 
use the values of the parameters m and g which provide 
the best fits to the lattice results in the Landau gauge in 
Ref. [41], that is, m = 0.54 GeV and g = 4.9 in the zero- 
momentum scheme and m = 0.39 GeV and g = 3.7 in 
the infrared-safe scheme. Notice that the relevant expan¬ 
sion parameter is 35 f^/( 167 r^) < 1. We have no reason, a 
priori^ to exclude a dependence of the parameters m and 
^ at a certain scale with the value of ^ at the same scale. 
Such dependences would have to be inferred from fits to 
lattice data. In the absence of such data, we assume fixed 
values of m and g adjusted from lattice data at ^ = 0 . 


Defining the renormalized coupling constant from the 
vertex Tacc{Pi ^—p)^ this relation holds to the finite parts 
as well, which completely fix the factor Zg. Unlike in the 
Landau gauge, where the relation (76) holds to all orders 
[56] . in the present case we have only checked that this 
is so at one-loop order. In particular, we do not know 
whether Eq. (76) is compatible with the renormalizabil- 
ity of the theory in general. But, as in the case of the 
second relation ( [T^ , we can always impose the prescrip¬ 
tion Eq. (76) for the finite parts only. As discussed in 
Ref. [24], the explicit calculation of the divergent part 
of Zg at one-loop order shows that the latter is both n 
and ^ independent and thus coincides with the standard 
result 


= -Pk, (77) 

yielding the usual universal one-loop beta function, as 
it should. More generally, we show in Appendix [P] that 
the present model reproduces the known one-loop UV 
behavior of Yang-Mills theories in the Landau gauge. 

The identity ( [T^ is valid only in the limit n ^ 0 with 
the definitions ( |45| ). It does not hold in the CE case 
(n = 1 ), where the contributions to the ghost-antighost- 
gluon vertex proportional to do not vanish. Eor the 
purpose of comparing the n ^ 0 and n = 1 results, 
we employ a renormalization scheme for n = 1 as close 
as possible to that used for n ^ 0. It is based on a 


A. Zero-momentum scheme 


Eigme[^ shows the ghost and transverse gluon propa¬ 
gators renormalized within the zero-momentum prescrip¬ 
tions ([5^-([57|) for different values of the gauge-fixing pa¬ 
rameter At small we recover the Landau gauge 
results of Ref. [41], where, in particular, the ghost prop¬ 
agator diverges at low momenta but the ghost dressing 
function p‘ ^G^ h{p) remains finite. The renormalization 
condition ( j^ imposes that for ^ 7 ^ 0 the ghost propa¬ 
gator itself remains finite in the infrared, where we see 
that it essentially behaves as a massive propagator. The 
strong ^ dependence at low momentum is governed by the 
second condition (53). As for the transverse gluon prop¬ 
agator, its values at p = 0 and p = p are constrained by 
om choice of renormalization conditions ([5^ and (54), 


which do not involve thus leaving little room for a 
possible ^ dependence. Indeed, in this scheme, the lat¬ 
ter only appears through loop effects and is thus small. 
Still, we observe that when ^ grows, the gluon propa¬ 
gator flattens near its p = 0 value. In fact, the lat¬ 
tice results in the Landau gauge show that the trans¬ 
verse gluon propagator is nonmonotonous at small p. Eor 


® The definition of the renormalization factor used here differs 

from that of Ref. ESI- 

® Results for SU(2) are qualitatively similar. 
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0.0 0.5 1.0 1.5 2.0 

p [GeV] 

FIG. 7: The ghost (top) and the transverse gluon (bottom) 
propagators as functions of momentum for various values of 
^ in the zero-momentum renormalization scheme with m = 
0.54 GeV and g = 4.9. 


^ = 0, we reproduce the one-loop behavior of Ref. lU: 

(p) ~ + ^^A^/(1927r^)p^ ln(p^//i^) + 0{p‘^). We 

observe in Fig. that this nonmonotorious behavior is 
more pronounced for increasing 

The above results show that the (renormalized) propa¬ 
gators are essentially governed by their tree-level expres¬ 
sions, respectively 

GT(P)= 2 A 2 Gr^(p) = ^^. (79) 

In order to emphasize the size and the shape of loop 
corrections, we plot in Fig. [^the loop contribution to the 
inverse propagators as a function of p^. We see that loop 
corrections are at the percent level for the ghost and at 
the ten percent level for the transverse gluon. 

B. Infrared-safe scheme 

The ghost and the transverse gluon propagators in the 
infrared-safe scheme are shown in Fig. Again, we 

check that we recover the results of Ref. [41] (in the 
appropriate renormalization scheme) for small values of 




{? [GeV^] 


FIG. 8: The loop contributions — (p^ + ^m^) (top) 

and G^^(p) — (p^ + m^) (bottom) to the ghost and the trans¬ 
verse gluon inverse propagators as functions of p^ in the zero- 
momentum scheme with m — 0.54 GeV and g — 4.9. 


In the present case, the value of the gluon propagator 
at vanishing momentum is not fixed and varies strongly 
with This can be traced to the fact that ^ influences 
the propagator at tree level through the second equation 
in Eq. ( [7^ . As in the previous case, we observe that the 
approach to the p = 0 value flattens as ^ is increased. 
We plot the loop contribution to the inverse propagators 
in Fig. These are again at the percent level for the 
ghost and at the ten percent level for the gluon. 


C. Comparison with the CF model 

We compare the results of the previous section to those 
of the CF model, obtained by setting the number of 
replica n = 1. As already mentioned, this allows us to 
pinpoint the peculiar effects of our treatment of the Gri¬ 
bov ambiguities. In particular, taking the limit n ^ 0 is a 
crucial step for our approach to correspond to a bona fide 
gauge-fixing procedure. We have already discussed the 
fact that an important difference with the case n = 1— 
already present at tree level—is the fact that the gluon 
propagator is exactly transverse in the case n ^ 0, in 
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0.0 0.5 1.0 1.5 2.0 

p [GeV] 

FIG. 9: The ghost (top) and the transverse gluon (bot¬ 
tom) propagator as a function of momentum for various val¬ 
ues of ^ in the infrared-safe renormalization scheme with 
m = 0.39 GeV and g = 3.7. 


contrast to the case n = 1. In the present section, we dis¬ 
cuss the differences between the two cases arising from 
the one-loop contributions to the ghost and the gluon 
propagators. 

The case of the zero-momentum scheme is illustrated 
in Fig. EH We observe that the ghost propagators for 
n ^ 0 and n = 1 are essentially the same while the 
transverse gluon propagators are qualitatively different. 
Contrarily to the former, the latter receives a direct con¬ 
tribution from a loop of superfield, proportional to n — 1, 
as discussed previously; see Eq. ( [^ . We see that the 
flattening near p = 0 described above in the case n ^ 0 
is absent in the CF model. For completeness, we show the 
momentum dependence of the longitudinal gluon propa¬ 
gator in the CF model in Fig. 

A similar comparison is made in Figs. and 14 for 
the infrared-safe scheme. The previous observations hold 
here too where the differences between n ^ 0 and n = 1 
are even more pronounced. For this reason, we show more 
values of ^ than in the previous case. Again we observe 
that the flattening of the gluon propagator near p = 0 
observed previously is absent in the CF model. Also we 
have checked that the nonrenormalization property (50) 




FIG. 10: The loop contributions Ggj^^(p) — (p^+^m^) (top) and 
C^^(p) — (p^ + m^) (bottom) to the ghost and the transverse 
gluon inverse propagators as functions of p^ in the infrared- 
safe scheme with m = 0.39 GeV and g = 3.7. We observe that 
the loop contributions to the ghost inverse propagator vanish 
at p = 0, as a consequence of the non renormalizati on t heorem 
discussed in the text below Eq. (75); see also Eq. (50). 


of the ghost vertex at zero momentum indeed holds in 
the case n ^ 0 but not for n = 1. Einally, the one-loop 
longitudinal gluon propagator of the CE model in the 
infrared-safe scheme is shown in Eig. 


VI. RENORMALIZATION-GROUP 
IMPROVEMENT 


We now take into account the RG running of the pa¬ 
rameters m^, p, and ^ and we investigate the RG tra¬ 
jectories in the two renormalization schemes discussed 
above. As usual, we define the beta functions [3^ for the 
parameters a = m?^g^ ^ as 


din Zoc 
d In /i Q ’ 

and the gluon anomalous dimension as 


Pa = 


da 


dln/i 


(80) 


din Za 
dinp Q 


( 81 ) 
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FIG. 11: The ghost (top) and transverse gluon (bottom) prop¬ 
agators computed in either the present gauge fixing (n ^ 0) 
or the CF model (n = 1) in the zero-momentum scheme with 
m = 0.54 GeV and g = 4.9, for ^ = 1 (left) and ^ = 4 (right). 


FIG. 13: The ghost propagator computed in either the present 
gauge hxing (n ^ 0) or the GF model (n = 1) in the infrared- 
safe scheme with m = 0.39 GeV and g = 3.7 for various values 
of 




FIG. 12: The longitudinal gluon propagator in the GF model 
(n = 1) as a function of momentum for various values of ^ in 
the zero-momentum scheme with m = 0.54 GeV and g = 4.9. 


FIG. 14: The transverse gluon propagator computed in either 
the present gauge fixing (n ^ 0) or the GF model (n = 1) in 
the infrared-safe scheme with m = 0.39 GeV and g — 3.7 for 
various values of 


where the subscript 0 means that the right-hand side is 
evaluated at fixed bare quantities. Note that the pre¬ 
scription (76) at one-loop order implies that the ghost 
anomalous dimension is given by 


din Zc 
dln/i Q 


'oTa H -• 

2 g 


(82) 


Similarly, the relation ( [7Q| ) which, as discussed above, is 
valid in the two renormalization schemes studied here, 
implies 


integrate numerically the flow equations ([8Q|)-([8^ with 
initial conditions at the scale /io = 1 GeV. We use as 
initial conditions the values of the mass and coupling 
parameters of the previous section, namely = 

0.54 GeV and g{jiio) = 4.9 for the zero-momentum scheme 
and m(/io) = 0.39 GeV and g{jiio) = 3.7 for the infrared- 
safe scheme, and we vary the gauge-fixing parameter 


7c 




(83) 


These RG functions are obtained from the expressions of 
the renormalization factors in each of the renormaliza¬ 
tion schemes described previously, that is from Eqs. ([^- 
(63) in the case of the zero-momentum scheme or from 


the corresponding ones in the infrared-safe scheme. We 


A. Zero-momentum scheme 


In this scheme, the first condition (53) implies the fur¬ 
ther relation 


Pr, 


1A= 2 > 


(84) 
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FIG. 15: The longitudinal gluon propagator in the CF model 
(n = 1) as a function of momentum for various values of ^ in 
the infrared-safe scheme. 


which, when combined with Eqs. (|82[) and (83), yields 


^ gjiA ( 

C(mo) fl(Mo) V / 


(85) 


The running of the parameters m, and obtained 
from the direct integration of the flow equations in the 
zero-momentum scheme, is shown in Fig. We find a 

Landau pole at a finite scale, where the parameters g and 
m diverge, for all values of ^(/io)- In the case ^ = 0, it 
was argued in Ref. [41] that this originates from the fact 
that, at small /i, the renormalization conditions (53) and 


(54), which imply that Gt(0) > Grig) conflict with the 
gluon propagator computed at one-loop order, which is 
an increasing function of p at small p. The Landau pole 
prevents one from following the RG flow down to deep 
infrared scales. The proposal of Ref. [41] to bypass this 
problem is to freeze by hand the RG flow under a certain 
arbitrary scale M, chosen above the Landau pole such 
that the latter is avoided. In the case C = 0? treated in 
Ref. [41], this could be justified by a detailed analysis of 
higher loop contributions in the infrared. We shall not 
attempt to generalize this discussion to the case C 7^ 0 
here. Still, we present for the sake of illustration, a modi¬ 
fied zero-momentum scheme where we freeze the RG flow 
by hand under a certain scale, in the spirit of Ref. [41]. In 
practice, we use an RG scale p = ^/p^ + M^, where p is 
the actual momentum scale at which we evaluate corre¬ 
lation functions. In what follows, we use M = 0.5 GeV, 
which is of the order of the typical value of the mass 
parameter m. The corresponding flow of the coupling 
constant as a function of p is shown in Fig. [^ 

We implement this modified RG scheme to compute 
the flow of all parameters from which we evaluate the 
RG-improved ghost and gluon propagators, shown in 
Fig. We see that, despite the relatively important 
flow of the parameters, the propagators are almost unaf¬ 
fected as compared to those of strict perturbation theory 
in this scheme; see Fig. [^ 



o'—■—■—■—■—^—■—■—■—■—■—■—■—■—■—■—■—■—■—■—' 

0.0 0.5 1.0 1.5 2.0 



A/ [GeV] 



FIG. 16: The running of the parameters g{p), Tn{p), and ^(p) 
in the zero-momentum scheme, for various values of ^(/xo) = 
We observe a Landau pole at a finite scale, where g and m 
diverge, and the RG flow stops. 


Although the above procedure allows one to avoid the 
Landau pole, the need to introduce the arbitrary scale M 
is somewhat uncomfortable. Moreover, despite the fact 
that most excitations in the present theory are massive, 
there remain massless degrees of freedom in the replicated 
superfleld sector (which is only absent for the CF model, 
n = 1) and it is not clear that the freezing of the RG flow, 
implemented by hand in the above procedure, is justified. 
It is the purpose of the infrared-safe scheme to provide a 
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0.0 0.5 1.0 1.5 2.0 


p [GeV] 



FIG. 17: The parameters and ^(/x) as functions of 

the momentum p, with ji = and M = 0.5 GeV in 

the zero-momentum scheme, for various values of ^(/xo) = 
The RG flow is effectively frozen for p < M, which avoids the 
Landau pole. 



0.0 0.5 1.0 1.5 2.0 

p [GeV] 

FIG. 18: RG-improved ghost (top) and transverse gluon (bot¬ 
tom) propagators as functions of momentum in the zero- 
momentum prescription scheme with /x = y/p‘^ for var¬ 

ious values of ^(/xo) = 


avoids the conflicting prescriptions at small /x discussed 
above. As a consequence, we find that the infrared- 
safe RG flow can be integrated down to arbitrarily small 
scales /X, depending on the choice of initial conditions'^, 
as was first pointed out in Refs. @01111] in the case C = 0- 
In this scheme, the relation Za/Z^ = 1 dis¬ 

cussed above implies that 


7A = - 




from which we get, using Eq. 


82h and (83), 


an) _ gHu) /m(Mo)y 

^(Mo) V rnin) ) 


( 86 ) 


(87) 


more systematic procedure which avoids the Landau pole 
[4T] . as we now discuss. 


Figure shows the RG flows of the parameters m 
and ^ for various values of C(/^o)- We observe that both 


B. Infrared-safe scheme 


The relaxation of the first condition (53)—replaced by 


the second condition (73) in the infrared-safe scheme- 


10 Pqj. instance, it is clear that in the case m(/xo) = C(mo) = 0 
[which implies m(/x) = ^(/x) — 0 foi" o-H P o.nd thus corresponds to 
the standard (Faddeev-Popov) Landau gauge] one gets a Landau 
pole. 
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the coupling and the mass first increase for decreasing fi 
and are then attracted towards zero in the infrared. We 
also see that the maximal values of both parameters de¬ 
crease with increasing ^(yUo) and are therefore, maximal 
in the Landau gauge. In all cases we have considered, 
the coupling remains small enough for perturbation the¬ 
ory to be (qualitatively) meaningful: we recall that the 
relevant expansion parameter is 3^^/(167r^). Finally, we 
observe that the gauge-fixing parameter is first attracted 
towards zero as /i decreases but eventually diverges in 
the limit /i ^ 0. In particular, we find that the Landau 
gauge fixed point = 0) is unstable in the infrared. 



IJ [GeV] 



0.0 0.5 1.0 1.5 2.0 

IJ [GeV] 



IJ [GeV] 


FIG. 20: Running of the ghost mass parameter mgh(/i) = 
in the infrared-safe scheme, for various values of 

^(Mo) = 


shown in Fig. Using Eq. ( [S^ , we have 


^gh(M) ^ ^(/i) 

rrighil^o) rn{ij) ‘ 


(89) 


We observe that mgh(/i) is attracted towards a nontrivial 
fixed point in the infrared, which shows that g{ii) ^ m{g) 
and ^(/i) ^ l/m^(/i) ^ l/g‘^{g) when g ^ 0. 

An important remark to be made here is that the RG 
flows of the independent parameters g and m never freeze 
out despite the fact that most degrees of freedom of the 
theory are massive. This is due to the fact that there 
remains, in fact, massless excitations in the superfleld 
sector, as can be seen from the tree-level propagators 
(25) and (27). This is a particular feature of the present 
theory and, more precisely, of the way Gribov copies are 
handled. As we shall see below the RG flow of the CF 
model, where all degrees of freedom are indeed massive, 
is qualitatively different. This shows that the freezing of 
the flow introduced by hand in the previous scheme is 
unrealistic. In fact, the nonfreezing flow observed in the 
present scheme leads to dramatic RG corrections on the 
gluon propagator, as we now discuss. 

We show in Fig. 


21 


the RG-improved ghost and (trans¬ 
verse) gluon propagators, where we used, for the RG 
scale, g = p. Using Eq. (83), we obtain, for the RG- 
improved ghost propagator. 


Ggh(p) = 


»^gh(p) 1 

’Wgh(Mo) ’ 


(90) 


FIG. 19: Running of the parameters g{g), m{g), and ^{g) in 
the infrared-safe scheme, for various values of ^(go) = 

It is interesting to check the running of the ghost mass 
parameter, 

( 88 ) 


where mgh(/i) is defined in Eq. (88). Since the latter 
saturates to a constant value in the infrared, as discussed 
above, we check that the value of the ghost propagator 
at zero momentum does not receive any RG correction: 


G'gh(O) = 


1 

”^gh(Po)’ 


mgh(M) =to(m)vU^ 


(91) 
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as expected from the nonrenormalization relation ( [7Q| ). 
This is easily checked in Fig. The fact that the ghost 
propagator is fixed to its tree-level value both at p = 
0 and p = /io also explains (together with the relative 
smallness of loop corrections in the ghost sector discussed 
previously) that the RG-improved ghost propagator is 
very similar to the one without RG improvement. 



0.0 0.5 1.0 1.5 2.0 

p [GeV] 


where we have defined the dimensionless mass parameter 
m(p) = m(/i)//i. The fact that the ghost mass parameter 
(IS^ reaches a plateau for sufficiently small p and not too 
small ^(po) implies that the p ^ 0 behavior of the gluon 
propagator is governed by the function fh? (p)/ [1+m^ (p)], 
a monotonously 



IJ [GeV] 


FIG. 22: Running of the dimensionless mass parameter 
m(p) = m(p)/p in the infrared-safe scheme, for various values 
of ^(po) = i 


increasing function of fh?{p). The running of the lat¬ 
ter is shown in Fig. where we observe that it is 
a monotonously decreasing function of p and that, for 
^(po) not too small, it reaches an infrared fixed point, 
which is approached linearly. We conclude that for ^(po) 
not too small, 


Gt(p ^ 0) 


^(Mo) 

”^gh(0) 1 + "1^(0) 


[1 + 0{p )], 


(94) 


where the linear term in p is negative. 


FIG. 21: RG-improved ghost (top) and transverse gluon (bot¬ 
tom) propagators as functions of momentum in the infrared- 
safe scheme with p = p, for various values of ^(po) = 


Things are different for the gluon propagator. In that 
case, the relation (86) implies that the RG-improved 


transverse propagator reads 


Grip) = 


^(Mo) 


1 


^{p) P^ + rTn?{p)' 


(92) 


We see in Fig. that the RG effects strongly modify 
the gluon propagator as compared to its strict one-loop 
expression; see Fig.|^ In particular, both its value at p = 
0 and the way the latter is approached are dramatically 
altered, even for small values of ^(po)- The flattening 
near p = 0, observed in strict perturbation theory, is 
turned into a linear behavior. This can be understood as 


follows. We first rewrite Eq. (92) as 


Grip) — 


^(Po) m^jp) 

ml^^ip) 1 + rh?(p) ’ 


(93) 


C. Comparison with the CF model 

We compare our results for the RG flow of the present 
theory (n ^ 0) with that of the CF model (n = 1)^^. 
We only consider the infrared-safe scheme, for which no 
freezing of the flow in the infrared is put in by hand. The 
infrared safety of this scheme is also verified in the CF 
model, as can be seen in Fig.|^ Clearly, the RG flow of 
the CF model is dramatically different from that of the 
case n ^ 0. In particular, we observe that the RG flow 
freezes below a certain scale. This is to be expected since 
the CF model only contains massive degrees of freedom 
which decouple in the deep infrared. The only exception 
is the case ^(po) = 0, where the ghost becomes massless. 
In this case, the CF model is equivalent to the case n ^ 0, 


11 We recall that, in the CF model, our prescription for the renor¬ 
malized coupling constant is given by Eq. ( |78| >. 
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fj [GeV] 



FIG. 23: The RG flow of the parameters m(/x), and ^(/x) 

in the GF model (n = 1) with the infrared-safe scheme. 


as already discussed, and we see that the RG flow does 
not freeze in the infrared. As a consequence of the RG 
freezing for ^(/xq) 7 ^ 0 , the coupling and mass parameters 
g and m do not vanish in the infrared and the gauge¬ 
fixing parameter ^ does not diverge. All the parameters 
reach constant values. 

The RG-improved ghost and gluon propagators for 
the CF model are shown in Fig. in the infrared-safe 
scheme. We see little change as compared to the results 
from strict perturbation theory (Figs. [T^[T^, for ^ not 



0.0 0.5 1.0 1.5 2.0 


p [GeV] 

FIG. 24: The RG-improved ghost and gluon propagators in 
the GF model (n = 1) with the infrared-safe scheme. 


too small, despite the relatively important change of the 
running parameters in the range of momenta considered 
here. RG corrections seem more important for small 
However, the RG-improved results in the CF model are 
dramatically different from those of the theory consid¬ 
ered here (n ^ 0), in particular for the gluon sector. 
Altogether these qualitative differences illustrate the im¬ 
portant role played by the replicated superfield sector of 
the theory and, in turn, by the Gribov copies, in partic¬ 
ular in the infrared. 
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VII. CONCLUSIONS 

We have investigated Yang-Mills propagators in the 
class of nonlinear covariant gauges recently proposed in 
Ref. [24]. These are defined from the minimization of 
a functional [see Eq. (§] that generalizes the one used 
in the case of the Landau gauge. In particular, it is ex¬ 
pected that minimization algorithms routinely employed 
in the Landau gauge can be applied to the present pro¬ 
posal. This would open the way to lattice simulations in 
covariant gauges other than the Landau gauge. 

The continuum formulation of these gauges relies on 
a generalization of the averaging procedure over Gribov 
copies proposed in Ref. m- The standard Faddeev- 
Popov construction—which neglects the Gribov issue— 
produces the GFDJ action. Instead, our treatment of 
Gribov copies results in the GF action—a massive gener¬ 
alization of the GFDJ action—coupled to a replicated set 
of n —1 supersymmetric nonlinear sigma models. The dy¬ 
namics of the replica sector is such that the theory indeed 
corresponds to a gauge-fixing procedure in the physically 
relevant limit n ^ 0. Our gauge-fixed action is pertur- 
batively renormalizable in four dimensions, which allows 
for perturbative calculations. In the Landau gauge limit, 
this procedure allows one to reproduce, with a one-loop 
calculation, the lattice results obtained in the minimal 
Landau gauge, where one selects a unique Gribov copy 
in the first Gribov region [Hiiiaiii]. This suggests that 
for a certain range of values of the averaging parameter 
Po, all copies in the first region are equiprobable and av¬ 
eraging over them is essentially equivalent to picking up 
a unique one. If this property was to hold away from 
the Landau gauge limit as well—which can only be as¬ 
sessed by actual lattice calculations— we expect that the 
present results provide reasonable predictions for possible 
lattice calculations in the gauge proposed here. 

A key point of the continuum formulation concerns 
the replica technique and the interplay between the limit 
n ^ 0 and renormalization. Indeed, part of the n de¬ 
pendence of the theory can be absorbed in the definition 
of the renormalized parameters [24] |40] . In the present 
work, we have employed a minimal scheme which, first, 
allows one to reproduce the results of Ref. m in the Lan¬ 
dau gauge (^ = 0) and, second, has a smooth ^ > 0 limit. 
We have studied two renormalization schemes and have 
presented results for the ghost and gluon propagators at 
one-loop order with and without RG improvement. Fi¬ 
nally, we have compared our results to those of the GF 
model, obtained by simply setting the number of replicas 
n = 1. This allows one to pinpoint the peculiar effects 
of the superfield sector of our theory, which is related to 
our particular treatment of the Gribov ambiguities. 

The first important aspect of the present treatment of 
Gribov copies is the fact that the basic fields of the the¬ 
ory acquire effective masses, related to the gauge-fixing 
parameters /3o and ^o- In contrast to the Landau gauge 
(Co = 0) case, not only the transverse gluons, but also the 
FP ghosts are massive. A striking difference between the 


GF model (n = 1) and the gauge-fixed theory (n ^ 0) 
is the fact that, in the latter case, the gluon propagator 
remains transverse in momentum space even away from 
the Landau gauge. This is already visible at tree level. 
At one-loop order in a strict perturbative expansion— 
i.e., without RG improvement—we also observe impor¬ 
tant differences between the n = 1 and the n ^ 0 cases, 
mainly in the transverse gluon propagator, which receives 
direct contributions from the replica sector. It is impor¬ 
tant to note that these contain massless degrees of free¬ 
dom, which lead to nonanalyticities at small momentum. 
These are absent in the case n = 1. 

The role of massless excitations in the gauge-fixed the¬ 
ory is further illustrated by implementing RG improve¬ 
ment. We have devised an infrared-safe renormalization 
scheme, generalizing the one put forward in Ref. [41] 
for the Landau gauge. We observed that the one-loop 
flow has no Landau pole and that both the coupling and 
the mass parameters are attracted towards the Gaussian 
fixed point, whereas the Landau gauge fixed point ap¬ 
pears unstable, both in the UV and in the infrared. This 
can be seen as a direct effect of the presence of massless 
modes in the theory. In the Landau gauge limit (C = 0), 
these are the ghost fields whereas for C 7 ^ 0, these come 
from the replica sector. Still, we observe that the C ^ 0 
limit is smooth. This results in strong RG effects on the 
propagators, mainly the transverse gluon, down to the 
deep infrared. These results of the gauge-fixed theory 
are to be compared to the corresponding ones in the GF 
model, where all degrees of freedom are massive and the 
RG flow freezes in the infrared, resulting in quantitative 
but not qualitative changes as compared to the strict 
perturbative results. 

The present work provides an explicit example where 
one can explicitly work out the effect of the Gribov copies 
in a semianalytical continuum calculation, in a class of 
gauges that may be amenable to lattice calculations. In 
the case of the Landau gauge, our treatment of Gribov 
copies already allowed us to capture some nontrivial in¬ 
frared physics in perturbation theory, whereas standard 
techniques based on the FP quantization require the use 
of nonperturbative tools. However, apart from provid¬ 
ing an effective bare gluon mass, the replica sector com¬ 
pletely decouples in the Landau gauge. This is not so 
in the present case and one can thus study explicitly the 
role of the replica sector—and thus, indirectly, of the 
Gribov copies—in the dynamics of the theory in the in¬ 
frared. We believe the present results are of interest for 
the general question of possible gauge (in)dependences 
of the infrared sector of the Yang-Mills propagators. We 
hope that the present work will stimulate lattice studies, 
e.g., along the lines of Refs. [T7ll2Q] . 
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Appendix A: Complete propagators 

We detail the inversion of the quadratic part of the ac¬ 
tion in the sector (A,ih,A/c). The most general form 
of the various two-point vertex functions is given in 


Eqs. (30)-(33). These are constrained by the spacetime 
symmetries, the replica permutation symmetry and the 
isometries of the Grassmann subspaces for each replica. 
These can be grouped together in the matrix 


T^^^ = I ihA ^ihih 0 

\ 0 (5fci + r2n0j, + ((5fe; - i)r3^ 


(Al) 


where the scalar functions Fr, and only depend on and where is the Laplace operator on the 

curved Grassmann space, defined as 1211150 ] 


U0_ = = 2l3o{ed0 + Ode) + 2(1 - Me)ded-e. 

V9{0 


(A 2 ) 


In particular, it satisfies the identity 0^6(0, O') = —2 + ‘^l^oS{0,0'). 

Note that, since there are n — 1 replicas, the matrix representation only makes sense for n > 0 and the limit 
n ^ 0 must be done after the inversion. Using the symmetries of the problem , the most general form of the inverse 
matrix reads 


+ ^hPpy iPp\hA -iPp^A 

1 _ I A5 

ipv^A A 5 5im [Ai(i {0_l,0^) + {0i,0_m)\ + (1 ~ ^Im) Asy 


(A3) 


where the unknown scalar functions A^, A/,, and Ai ^...^5 only depend on The inversion is defined by 


0 


r( 2 )x = (T 1 0 

^ ^ '00 6kmS{0^^0J^ 


(A4) 


where the product x involves a sum over Lorentz and replica indices and an integral over the Grassmann variable 0i . 
The calculation is straightforward. We get, for the gluon propagator. 


At = and Ar = 


^ihih (ri P {n — 2)^or3) 


{rLTihih - P^r^ihA) (ri + - 2 )/ 3 or 3 ) - (n - ' 


The other components in the (A, ih) sector are obtained from these as 


AihA = 


- ihA 


^ihih 

Finally, the components in the superfield sector read 

A4 = - ^ 


and Aihih = 


1 P p‘^rihAAihA 


- ihih 


As = -p 
As = 


(n - l)/3op2F4 

A 4 




^ihAiJ^L^ihih P 

^ihA 


2 r ihA 


^ihih 

1 


Fi + (n - 2)/3 oF3 
Al = -h /3 oA3 


A 2 = 


Ti - F3/3 o 

1 


P^F4A4 - 


Al 


(Fi - /loFs), 


2/3o(Fi + 2^or2) 2^0 ■ 


(A5) 

(A 6 ) 

(A7) 

(A 8 ) 

(A9) 

(AlO) 

(All) 


The propagators are obtained by taking the limit n —> 0. by 
For instance, the longitudinal gluon propagator is given 


Gl{p) = lim Ai(p), 


(A 12 ) 
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which is Eq. (36). 


and by performing a Legendre transform with respect to 
the sources 77 ^, 77 ", and M". It reads 


Appendix B: Tree level 


We give here the tree-level expressions of the bare two- 
point vertex functions. In the transverse gluon and the 
ghost sectors, 


rr(p) =p^ + npo 

(Bl) 

rcc(p) =p^ + ^ 06 ; 

(B2) 

In the longitudinal gluon and h sector. 

rL(p) = npo 

(B3) 

^ihAip) = 1 

(B4) 

^ihih{p) = -^ 0 ; 

(B5) 

In the superfleld sector 

ri(p) =p^ 

(B 6 ) 

r 2 (p) = Co /2 

(B7) 

r3(p) = 0 

(B 8 ) 

r4(p) = 1 . 

(B9) 


We check that the Slavnov-Taylor identities of Ap¬ 
pendix are satisfied at tree level. We also check that 
the right-hand side of Eq. (36) yields (for arbitrary n) 


Gl{p) = 


Co 


+ /3oCo' 


(BIO) 


Appendix C: Slavnov-Taylor identities in the CF 
model 

As already discussed, the CE model can be obtained 
from the theory considered here by setting the number 
of replicas n = 1. Its action is given by 

S[A,c,c,h]=SYu[A] + ScF[A,c,c,h] (Cl) 

with S'ym and 5 'cf given in Eqs. 0 and ( p^ respectively. 
This model possesses a (non-nilpotent) BRST symmetry, 
whose action on the fields is 

sAl = D^c\ sc“ = -|r'’^cV (C2) 

and 

sc^ = ih^, sih^ = Poc^. (C3) 

The Zinn-Justin equation corresponding to this symme¬ 
try is obtained as usual, i.e., by introducing external 
sources for the fields and for all the independent BRST 
variations, S ^ S — Si with 

(C4) 


f ( 5T 5T (5r (5r „„ (5r „ „ (5r 1 „ 

L 1 6K^ SAf^ (5X“ (5c“ ^ (5c“ Sih°- j ~ 

(C5) 

Taking two derivatives of this equation with respect ei¬ 
ther to ih^ and or to and and setting the sources 
to zero, one obtains the following symmetry identities for 
the two-point vertex functions in momentum space, as 
defined in Eq. ( [^ : 

^cK^(.P)^A^ih(.p) - rcc(p) - Po^ihihip) = 0 (C 6 ) 

and 

Tc/?, {p)^A,A. ip) - P^A^ihA. (p) = 0. (C7) 

Eliminating E^^^ (p), we obtain, finally, 

p' [TihA{p)? - TL{p)Tihih{p) = LkMd, (C 8 ) 

Po 


which is used in obtaining Eq. (38). 


Appendix D: RG flow in the ultraviolet 


The flow of the present theory in the ultraviolet is eas¬ 
ily obtained from the divergent parts of the renormaliza¬ 
tion factors; Eqs. (64)-(69). Eor instance, the /3 functions 
for the coupling and the gluon mass parameters read, in 
the limit p ^ m, rugh, and taking n ^ 0 , 


yp 

9 


11 g'^N 
y 871 ^ ’ 


and 


3UV 




35 g'^N 
12 8^"' 


(Dl) 


The UV running of the coupling is the same as in the 
Landau gauge, that is 


1 _ 1 IIN 

g^(fi) g^ip*) 247r2 V/x*/ ’ 


(D2) 


where /i and p* are UV scales. The flow of m is easily in¬ 
tegrated. Using Eqs. (87) and ( 88 ), we get the following 
relations 


min) ^ / g{p) 

rrighip) ^ f gjp) W 
TOgh(p*) \g{p*)) 

C{p) ^ / g{p) Y“ 
C(p*) \g{p*)) 


(D3) 

(D4) 

(D5) 


We observe that the Gaussian point p = 0 and m = 
rugh = 0 is an attractive UV fixed point, whereas C = 0 
is an unstable UV fixed point. 
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Using Eqs. (90) and (92), we recover the standard one- 


loop behaviors of the gluon and ghost propagators at 
large momentum 


where we have used the fact that the contribution 
from the mass terms are negligible at large momentum: 


(P) ^ 


G^'^ip) (X 


1 




(D6) 

(D7) 
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